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ABSTRACT 



Aims. We report on the detection of vinyl cyanide (CHiCHCN), cyanomethyl radical (CH 2 CN), methylacetylene (CH 3 CCH) and 
thioformaldehyde (H2CS) in the C-rich star IRC +10216. These species, which are all known to exist in dark clouds, are detected for 
the first time in the circumstellar envelope around an AGB star. 

Methods. The four molecules have been detected trough pure rotational transitions in the course of a A 3 mm line survey carried 
out with the IRAM 30-m telescope. The molecular column densities are derived by constructing rotational temperature diagrams. A 
detailed chemical model of the circumstellar envelope is used to analyze the formation of these molecular species. 
Results. We have found column densities in the range 5 x 10 12 - 2 x 10 13 cirT 2 , which translates to abundances relative to H2 of 
several 10~ 9 . The chemical model is reasonably successful in explaining the derived abundances through gas phase synthesis in the 
cold outer envelope. We also find that some of these molecules, CH 2 CHCN and CH 2 CN, are most probably excited trough infrared 
pumping to excited vibrational states. 

Conclusions. The detection of these species stresses the similarity between the molecular content of cold dark clouds and C-rich 
circumstellar envelopes. However, some differences in the chemistry are indicated by the fact that in IRC +10216 partially saturated 
carbon chains are present at a lower level than those which are highly unsaturated, while in TMC-1 both types of species have 
comparable abundances. 

Key words, astrochemistry - circumstellar matter - stars: AGB and post-AGB - stars: carbon - stars: individual: IRC +10216 



1. Introduction 

The chemical complexity of circumstellar envelopes (CSEs) 
around carbon-rich AGB stars is illustrated by the well known 
object IRC +10216, where more than 60 different molecular 
species have been identified to date. Due to the carbon-rich na- 
ture of the gas, the vast majority are oxygen-deprived hydro- 
carbons with a very concrete structure consisting of a linear 
and highly unsaturated backbone of carbon atoms. Among these 
species are cyanopolyynes HC2, !+ iN, polyacetylenic radicals 
C„H, carbenes H2G,, allenic radicals HC2«N as well as sulphur 
and s ilicon-bearing carbon chains C„S, SiC„ dCernicharo et al.L 
2000). The synthesis of these molecules involves the photodisso- 
ciation and photoionisation by interstellar UV photons of parent 
species out-flowing from the star and subsequent neutral-neutral 
and i on-molecule gas phase reactions in the co l d outer envelope 
(e^LafWeTaflll98a iNeiad & Millanfl98H IMiilar & Herbsll 



however a sizable fraction of partially saturated species , 
such as methyl p olyyne s CH3C2„H (n=l,2 llrvine et aTT Fl98 1 : 
IWalmslev et al.l Il984l). me th ylcyano p olyynes CH 3 Cja+iN 



yi' 

(n=0. 1 .2 iMatthews & Searsi Il983at iBroten et al.l 11984 



Snyd er et al. 1 120061) . cyanomethyl r adical CH2CN ( Irvine et al 



1988), vinyl cyanide CHIC HCN dMatthews & SearsL 11983b) 



1994; IMiilar etldl 120001) . Thus, reactive molecules are dis- 
tributed in circums tellar shells, as co nfirmed by interferometric 
observations (e.g. lGuelin et al.fl993l) . 

This type of chemistry resembles that occurring in cold 
dense clouds, such as TMC-1, which are also rich in highly 
unsaturated carbon chain molecules. The unsaturated character 
is typically the result of low temperature non-equilibrium 
chemistry and reflects the low reactivity of H2 with most 
hydrocarbons at low temperatures. Dark clouds contains 



cyanoallene CH 2 CCHCN (lLovas et all 120061) ; and highly 
saturated hydrocarbon s such as propylene CH2CHCH3 
dMarcelino et al.Ll2007l) . 

The question of why those partially saturated species are not 
detected in IRC +10216 is a critical one for the understand- 
ing of low temperature non-equilibrium chemistry. We have 
thus embarked on a deep search for partially saturated organic 
molecules in IRC +10216. In this paper we present the detec- 
tion of CH 2 CHCN, CH 2 CN and CH 3 CCH, which are observed 
for the first time in a CSE around an AGB star. The related 
species CH3CN has b een already identified in this source by 
iJohansson et al.l dl984l) . Finally, we also report on the detection 
of thioformaldehyde in IRC +10216. 

2. Spectroscopy and observations 

The four molecules on focus in this paper (CH2CHCN, CH2CN, 
CH3CCH and H2CS) were i dentified in space for the first 



time towards Sag i ttarius B2 dGardner & Winnewisset 



Send offprint requests to: M. Agundez 

* Based on observations carried out with the IRAM 30-meter tele- 
scope. IRAM is supported by INSU/CNRS (France), MPG (Germany) 
and IGN (Spain). 



llrvine et all 119881: ISnvder & Buhll [19731; ISinclair etaL 
by observation of rotational transitions. The microwave spec- 
trum of all these species has been extensively studied in the 
laboratory so that their spectroscopic properties are accurately 
known; see a compilation in the JPL and Cologne databases 
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Fig. 1. Rotational lines of (a) CH 2 CHCN, (b) CH 2 CN, (c) CH 3 CCH, (d) CH 3 CN and (e) H 2 CS observed toward IRC +10216 with 
the IRAM 30-m telescope at a spectral resolution of 1 MHz. The frequency scale is computed for an LSR source velocity of -26.5 
km s -1 . For comparison purposes, the thick grey panels show spectra corresponding to CRL 618 (from IPardo et al J 1200 7). "U" 
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Table 1. Line parameters of CH 2 CHCN, CH 2 CN, CH 3 CCH and H 2 CS in IRC +10216 
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Observed 


Calculated 
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Symmetry 


Frequency 
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Strength 


F 

'-'up 










(MHz) 


(MHz) 


(km s _1 ) 




(K) 


(Kkms" 1 ) 




CH 2 CHCN 


9o.9"8o,8 




84946.7(3) 


84946.000 


15.5(10) 


9.00 


20.4 


0.11(2) 


0.82 


"1,9-81,8 




83207.5(8) 


83207.505 


14.8(10) 


8.89 


22.1 


0.14(2) 


0.82 


92.7-82,6 




85715.4(5) 


85715.424 


14.5 d 


8.56 


29.2 


0.09(3) c 


0.82 


93,7-83,6 




85427.9(10) 


85426.920 


U.5 d 


8.00 


40.0 


0.09(3) r 


0.82 


93,6-83,5 




85435.6(10) 


85434.526 


14.5 rf 


8.00 


40.0 


0.07(2) c 


0.82 


10o,10"9o,9 




94276.3(10) 


94276.634 


U.5 d 
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25.0 


0.15(3) 


0.81 


10i,io-9i, 9 
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9.90 


26.6 
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0.81 


10 2 ,9-9 2 , 8 
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0.81 


10 2 ^-9 2 ,7 
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14.5 d 


9.60 


33.8 


0.14(3) 


0.81 


103,8-93,7 




94926.2(15) 


94928.606 


14.5 rf 


9.10 


44.5 


0.09(3) 


0.81 


10 3 ,7-9 3 ,6 




94942.0(10) 


94941.630 


14.5 rf 


9.10 


44.5 


0.08(2) 


0.81 


10 4 ,7-9 4 ,6 




94913.4(10) 


94913.115 


14.5 rf 


8.40 


59.7 


0.11(3) r 


0.81 


10 4 ,6-9 4 ,5 






94913.226 




8.40 


59.7 


C 


0.81 


Hi.ii-lOi.io 




101637.1(8) 


101637.231 


14.5 rf 


10.9 


31.5 


0.07(3) 


0.80 


112,10-102,9 




104211.5(10) 


104212.646 


14.5 rf 


10.6 


38.7 


0.08(3) 


0.79 


CH 2 CN 


4 , 4 -3o,3 J=9/2-7/2 


O 


80480.5(10) 


80480.384 


14.5" 


4.44 


9.7 


0.23(6) c 


0.82 


4o, 4 -3o,3 J=7/2-5/2 


O 


80489.6(10) 


80490.261 


14.5 rf 


3.43 
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O 
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P 
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14.5 d 
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0.82 


4i i3 -3i, 2 1=9/2-7/2 


P 


81232.4(5) 


81232.654 
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0.09(2) 


0.82 
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P 


99689.5(8) 


99689.833 


14.5 rf 


5.23 


13.4 


0.14(4) 


0.80 


5i,4-4i, 3 1=11/2-9/2 


P 


101531.6(10) 


101532.055 


14.5 d 
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0.22(6) 
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C 
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17.2 


0.22(5) c 
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J K =6 l -5 l 






102546.023 




5.83 


24.4 


C 


0.80 


H 2 CS 
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2.67 


8.1 
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Note.- The numbers in parentheses are errors in units of the last digits; the lcr error in J~T^dv is derived from least-squares 
fits and does not include the calibration uncertainty of 10 %. " v„ p stands for expansion velocity and is computed as half 
the full linewidth at zero level. b rj b is B cff /F cff , i.e. the ratio of T* A to T mb . A superscript "c" denotes a partly blended line 
while "d" indicates that the linewidth parameter \ exp has been fixe d to 14.5 km s -1 . T he calculated frequencies of CH 2 CN 
correspond to the strongest hyperfine component, as measured bv lOzeki et alj d2004). CH 2 CN and H 2 CS lines have ortho 
(O) or para (P) symmetry; E„,, is the energy of the upper level of the transition above the O or P ground state. 



for m olecular spectroscopy (Pic kett et all 119981 : iMuller et"aTL 
120051) . 

CH 2 CHCN is a planar asymmetric rotor with a- and b- 
type allowed tr ansitions (fi a = 3.815 D and = 0.894 D 
jStolze & Suited Il985l) . All the transitions observed in IRC 
+ 10216 belong to the stronger a-type. 

The radical CH2CN has two interchangeable hydrogen nu- 
clei with non zero spin which result in two distinct groups of ro- 
tational levels: ortho (parallel spins) and para (antiparallel spins) 
with relative statistical weights 3:1, between which both radia- 
tive and collisional transitions are highly forbidden. The elec- 
tronic ground state is 2 B\, thus, the quantum number K a is even 
for ortho levels and odd for para levels. The para ground state 
(li,i) lyes 14.15 K above the ortho ground state (Oo,o)- This as- 
signment is reversed compared to more common species with a 



A] electronic ground state (H 2 CO, H 2 CS, ...), where the quan- 
tum number K a is odd for ortho levels and even for para levels, 
with the rotational ground state having para symmetry. It s dipol e 
moment is relatively large: jj a — 3.5 D (Oze ki et ail [2004). 
The unpaired electron causes spin-rotation coupling which splits 
each transition /V , ,-N „ „ into two components. Further hy- 

K a ,K C K a ,K C 

perfine coupling of rotation with the non-zero spin of 'H and 14 N 
nuclei produces a myriad of components wh ose frequencies ha ve 
been precisely measured in the laboratory dOzeki et alll2004l) . 

Propyne (CH3CCH) is a prolate symmetric top molecule 
whose rotational levels are divided into two different species, 
A-type and E-type, not connected radiatively. Its dipol e moment 
is relatively small: fi a = 0.784 D dBurrell et all Il980h. The ro- 
tational spectrum was first measured by iTrambarulo & G ordv 
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(1950b and is now known with a hig h accuracy for its ground 
and some excited vibrational states ( Miille /et all l2002h . 

H2CS is an asymmetric rotor which, analogously to CH2CN, 
has two interchangeable hydrogen nuclei so that its rotational 
levels are grouped into ortho (K a odd) and para [K a even), with 
statistical weights 3:1. The ortho ground state ( 1 1,1) lyes 14.9 K 
above th e para ground s t ate (O q.q). Its dipole moment was mea- 
sured bv lFabricant et all d 19771) to be /i a = 1.647 D. 

The observations shown in Fig.[TJwere made with the IRAM 
30-m telescope during several sessions from 1990 to 2006, most 
of them after 2002 in the context of a A 3 mm line survey of 
IRC +10216 (Cernicharo et al. in preparation). SIS receivers op- 
erating at 3 and 2 mm were tuned in single sideband mode, with 
typical image rejections larger than 20 dB at 3 mm and around 
15 dB at 2 mm. We express intensities in terms of T^, the an- 
tenna temperature corrected for atmospheric absorption and for 
antenna ohmic and spillover losses. The uncertainty in T* due to 
calibration is estimated to be around 10 %. Data were taken in 
the standard wobbler switching mode with a beam throw of 4' . 
Pointing and focusing were checked by observing nearby planets 
and the quasar OJ 287. The back end used was a filterbank with 
a bandwidth of 512 MHz and a spectral resolution of 1.0 MHz. 
The system temperature was 100-150 K at 3 mm and 200 K at 
2 mm. On source integration times ranged from 2 to 20 hours, 
resulting in rms noise levels of 1-3 mK in T^. 

3. Results 

Fig. [TJshows the observed spectra of IRC +10216 correspond- 
ing to the CH 2 CHCN, CH 2 CN, CH3CCH and H 2 CS lines. The 
observational line parameters, obtained by using the SHELL 
fitting routine of the CLASS packagfl are given in Table [TJ 
Linewidths have been fixed in many cases due to limited signal- 
to-noise (S/N) ratio or due to blending with other lines. The 
measured widths are consistent with an expansion velocity v exp 
around 14.5 km s _1 , in agreement with most of the molecu- 
lar lines arising from th e expanding envelope of IRC +10216 
(ICernicharo et all 12000). There are no missing lines of CH2CN, 
CH3CCH or H2CS in the 3 mm atmospheric window. However, 
some lines of CH2CHCN having similar strengths and upper 
level energies to those reported in Table [TJ are not listed due to a 
complete blending with a stronger line or due to a low sensitivity 
of the spectra. For example the 1 lo,n-10o,io transition at 103575 
MHz is blended with a strong component of the 7=21/2-19/2 
2 ni /2 v 7 =l doublet of C 4 H. 

The line profiles can give us information on the spatial distri- 
bution of the molecules. In a spherically expanding envelope and 
for lines with a low optical depth, a double-peaked shape indi- 
cates that the emitting species has a distribution with an angular 
extent comparable or larger than the telescope beam (HPBW = 
21 "-31" for IRAM 30-m at A 3 mm), while a flat-topped pro- 
file indicates that the emitting region is not spatially resolved 
by the telescope beam. The lines of CH 3 CN in IRC +10216 
have a double-peaked character (see Fig. [TJ, which indicates 
that this molecule has an extended distribution. In the case of 
the four species with which we are concerned, the low S/N ra- 
tio of most of the observed lines makes difficult to distinguish 
whether line profiles are predominantly double-peaked or flat- 
topped for a given species. Therefore, any conclusion about the 
spatial distribution of the molecules derived from the line pro- 
files must be taken with caution. We will, nevertheless, consider 
that CH 2 CHCN, CH 2 CN, CH3CCH and H 2 CS have an extended 
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Fig. 2. Rotational temperature diagrams of CH2CHCN and 
CH 2 CN in IRC +10216. T ra , and N /0/ have been determined 
from a weighted least-squares fit to the observational points, and 
their uncertainties reflect the lcr errors in the fits to the line pro- 
files plus the 10 % estimated uncertainty in the calibration. 



distribution, based on chemical arguments (see §3.2 and Fig. [3), 
which essentially fills the telescope beam. Under such hypothe- 
sis, we have constructed rotational temperature diagrams to de- 
rive beam averaged column densities. 

For CH2CHCN (see Fig. [2]i we derive a total column density 
N for = (5.5 ± 1.5) x 10 12 cirT 2 and a rotational temperature of 
46 ± 16 K, which is within the range of rotational temperatures 
deriv ed for other shell-distribu ted molecules in IRC +10216: 20- 
50 K (ICernicharo et all 120001) . 

In the case of CH2CN, we have used all the observed lines - 
ortho and para- in a single rotational diagram (see Fig . [2J . Given 
that the lines are detected with a low S/N ratio, we do not aim 
at determining the O/P ratio although the available data is con- 
sistent with the statistical value 3:1. We derive a rotational tem- 
perature of 50 ± 13 K and a total column density of N, 0/ = (8.6 
± 1.4) X 10 12 cm" 2 , which is comparable to that of the related 
species CH3CN (see Tabled. 

The two CH3CCH lines detected are very likely a sum of the 
K — 0, 1 components, although the low S/N ratio do not allow to 
clearly distinguish them. Assuming a rotational temperature of 
30 K we tentatively derive a total column density of 1.8 x 10 13 
cm -2 . This value is comparable to the column density of CH3CN 
in IRC +10216, which is detected through lines about 25 times 
stronger than those of CH3CCH due to its much larger dipole 
moment (ficHj.CN - 3.925 D versus hchjCCH - 0.780 D). 

Both CH 3 CN and CH3CCH have b een also detected in the 
C-rich preplanetary nebul a CRL 618 (ICernicharo et al.l [200 1 ; 
iPardo & Cernicharol 120061: iPardo et al.l [20071: see also FiglTV 
which is in a evolutionary stage immediately following that of 
IRC +10216. However, the CH3CCH/CH3CN ratio is noticeably 
different in these two objects. In IRC +10216 both molecules 
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Table 2. Column densities of selected molecules in IRC +10216 



angular distance (") 
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Fig. 3. Abundances of CH 2 CHCN, CH 2 CN, CH 3 CCH and H 2 CS 
(solid lines) and related species (dotted and dashed lines) given 
by the chemical model, as a function of radius (bottom axis) and 
angular distance (top axis) for an assumed stellar distance of 150 
pc. 



are present with a similar abundance (CH3CCH/CH3CN = 
0.6). In contrast, in CRL 618 m ethylacetylene is m uch more 
abundant than methyl cyanide (Pardo et al., 2007), not only 
in the region corresponding to the first post-AGB ejections 
(CH3CCH/CH3CN = 6) but also in the warm and dense inner 
regions where UV photons efficiently drive a rich C-based pho- 
tochemistry (CH3CCH/CH3CN = 15). Such a stage has not yet 
been reached by IRC + 1 02 1 6. 

The four lines of H 2 CS observed in IRC +10216 have simi- 
lar E up values, thus it is rather difficult to constrain the rotational 
temperature. Assuming an statistical O/P ratio of 3: 1 and a rota- 
tional temperature of 30 K we derive a total column density N, 0/ 
= 1 x 10 13 cm~ 2 for thioformald ehyde, which is almost twice 
larger than that of formaldehyde (Ford et al. 2004; see Tabled. 



4. Chemical modelling 

In order to explain how the detected species are formed we have 
constructed a detailed chemical model of the outer envelope. 
The chemical network consists of 385 gas phase species linked 
by 6547 reactions , whose rate constants have bee n taken from 
ICernicharol d20 04). Agu ndez & Cernicharol (2006) and from the 
UMIST Database for Astrochemistrv dWoodall et aUl2007l) . The 
temperature and density radial profiles as well as other physi - 
cal assumptions are taken from Agun dez & Cernicharol ((2006). 
The resulting abundance radial profiles for CH 2 CHCN, CH 2 CN, 
CH3CCH, H 2 CS and related species are plotted in Fig. [3] The 
model predicts that these four molecules would display an ex- 
tended shell-type distribution with an angular radius of about 
20". The predicted and observed column densities are summa- 
rized in Table [2] 

Vinyl cyanide results from the reaction between CN and 
ethylene, the latter being found in the inner envelope with a n 
abundance relative to H 2 of 2 x 10~ 8 dGoldhaber et all 119871) . 
This reaction has been studied in the laboratory a nd found to be 
very rapi d at low tempera ture dSims et all IT993) yielding vinyl 
cyanide dChoi et al .1120041) : 
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Our model predicts a total column density N f0 , = 1.4 x 10 13 
cm" 2 , which reasonably agrees with the value derived from the 
observations^ The analogue of reaction (1) replacing C 2 FLt by 
C 2 H 2 is the main formation pathway of HC3N in IRC +10216. 
The much larger column density of HC3N compared to that of 
CH 2 CHCN (more than two orders of magnitude; see Table|2| re- 
flects the large excess of acetylene over ethylene in the inner en- 
velope of IRC +10216. In dark clouds, reaction ( 1 ) is also proba- 
bly th e main formation route to vinyl cyanide dHerbst & L eung, 
I1990I) . 

The dissociative recombination (DR) of the ion CH3CNI-F 
is the major pathway to form both CH 2 CN and CH3CN: 



CH + ^CH+^CHt^CH 3 CNH + 



>CH 2 CN 
>CH 3 CN 



(2) 



The branching ratios of the channels giving CH 2 CN and CH3CN 
are not known and are assumed to be equal. In our model, the 
major destruction process for both CH 2 CN and CH3CN is the 
photodissociation by interstellar UV photons. T he photodisso- 
ciation rate of CH3CN has been calculated by Roberg e~et alj 
d!991 ) but that of CH 2 CN is unknown and we assume it to be 
the same as that of CH3CN. Since the chemistry of these two 
molecules is relatively well constrained in IRC +10216, they are 
formed by the same reaction and their destruction is dominated 
by just one process, we may use the observed CH 2 CN/CH3CN 
ratio to estimate the branching ratios of the DR of the CH3CNFF 
ion. We find an 80 % for the channel (CH3CN + H) and 20 % for 
(CH 2 CN + H 2 or 2H). This estimate will be strongly affected if 
the photodissociation rate of CH 2 CN is very different from tha t 
of CH3CN but not if, as suggested by iHerbst & Leungl {1990) 
and lTurner et al.l dl990h . CH 2 CN reacts with atomic oxygen, the 
abundance of which is low in IRC +10216 at the radius where 
CH 2 CN is present. 

The related species CH3C3N is formed by the same se- 
quence of reactions that produces CH3CN, just replacing HCN 
by HC3N. Its predicted column density is lower than that of 
CH3CN by more than one order of magnitude. However, since 
our model predicts less CH3CN than observed, the CH3 ion is 



CN + C 2 H 4 



CH 2 CHCN + H 



(1) 



2 The radial column densities given by the model are multiplied by a 
factor 2 to get the total column densities across the envelope N fof . 
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probably underproduced, the column density of CH3C3N is most 
likely underestimated. If CH3C3N was present with a column 
density of several 10 11 cm -2 , the expected brightness tempera- 
tures would be a few mK and therefore it could be detectable. 
In fact, our IRAM 30-m line survey at A 3 mm of IRC +10216 
(Cernicharo et al. in preparation) shows two unidentified fea- 
tures with T* A -2-3 mK at 86756 MHz and 103280 MHz that 
could correspond to the J = 21-20 and J = 25-24 transitions of 
CH3C3N respectively. The spectra covering other 3 mm transi- 
tions of CH3C3N have not enough sensitivity to confirm or dis- 
card the presence of this species. 

The synthesis of CH3CCH involves various ion-molecule re- 
actions with the dissociative recombination of the ions C3H+ and 
C4H5 as the last step. The model underestimates the observed 
CH3CCH column density by one order of magnitude, probably 
due to uncertainties and/or incompleteness in the chemical net- 
work. The heavier chain CH3C4H is predicted to have a column 
density even higher than that of CH3CCH. The smaller rotational 
constant B and larger rotational partition function of CH3C4H 
makes it less favorable for being observed at millimeter wave- 
lengths. However its larger dipole moment (fJ.cH 3 c 4 H - 1-207 D 
versus ficH^ccH = 0.784 D) could result in line intensities simi- 
lar to those of CH3CCH. Our IRAM 30-m A 3 mm data shows 
an unidentified feature with T* A ~3 mK at 81427 MHz that could 
be assigned to the J = 20-19 transition of CH3C4H. As occurs 
with CH3C3N, the detection of CH3C4H remains tentative since 
our current spectra covering other 3 mm lines is not sensitive 
enough. 

La stly, as previously discussed by Agu ndez & Cernicharol 
(2006), H2CS is formed by the neutral-neutral reaction S + CH3 
and the dissociative recombination of H3CS + . The predicted col- 
umn density is N to , = 1.3 x 10 12 cm" 2 , which is a factor 7 lower 
than the value derived from observations. We note that a sig- 
nificant fraction of both H2CO and H2CS could be formed on 
grain surfaces by hydrogenation of CO and CS respectively. It is 
remarkable that in IRC +10216 thioformaldehyde is more abun- 
dant than formaldehyde, despite the cosmic abundance of oxy- 
gen being 50 times larger than that of sulphur. The carbon-rich 
character of the CSE makes oxygen-bearing species to have a 
very low abundance. 



5. Excitation mechanism: infrared pumping 

The chemical model predicts that the molecules formed in the 
outer envelope have their peak abundances between 3 x 10 16 cm 
and 6 x 10 16 cm, which corresponds to angular radii of 13"- 
26" for an assumed distance to IRC +10216 of 150 pc (see 
Fig. 0. Interferometric observations at millimete r wavelengths 
have located the molecular shell at about 1 5 ' ' dGuelin et all 
119931: lAudinos et all 1 19941 iLucas et all 1 1995b although the ex- 
act value depends on the molecule and transition mapped. At 
this distance the gas density is a few 10 4 cm" 3 and the g as ki- 
netic temperature is lower than 20 K ( Skin ner et all 1999). We 
may ask ourselves how the rotational levels are excited to result 
in rotational temperatures for CH 2 CHCN and CH 2 CN of 40-50 
K, well above the gas kinetic temperature. 

Collisions with H2 molecules do not seem to be respon- 
sible of such an excitation. For example, in the case of the 
cyanomethyl radical the collision coefficients y„i for CH2CN- 
H2 have not been measured or calculate d but should be simi- 
lar to those calculated for CH 3 CN-H 2 by iGreenl d!986l) . at least 
for the ^0=0 ladder (see lTurner et alj ( 1 19901) for more details). 
Thus, adopting a deexcitation collision coefficient y u \ = 10~ 10 
cm 3 s" 1 , typical for AJ = -1, AK = transitions of CH3CN, 



and an Einstein coefficient for spontaneous emission of A u \ = 5 
x 10" 5 s _1 , the resulting critical density is n crlt - A u i/y u i — 5 x 
10 5 cm" 3 . This value is larger than the gas density expected at 
a distance of 4-5 x 10 16 cm, so that the rotational levels are not 
thermalized and collisional excitation by itself should result in a 
rotational temperature lower than the gas kinetic temperature. In 
the cases of molecules such as CH2CN, with a relatively large 
dipole moment, i.e. whose rotational levels are hardly thermal- 
ized in the outer envelope, and which have rotational tempera- 
tures above the kinetic temperature of the gas, the excitation is 
most probably dominated by a radiative mechanism. 

Absorption of infrared photons and pumping into excited vi- 
brational levels followed by radiative decay to rotational levels 
in the ground vibrational state has been invoked many years 
ago as the main excitation mechanism fo r some circumstellar 
molecules in IRC +10216 dMorrisl [19751) . As a matter of fact, 
several molecules present in the outer envelope have been de- 
tected in excited vibrational states, e.g. C4H vj = 1, 2 (V7 = 
131 cm' 1 ; iGuelin et alJfl987T: lYamamoto et al]fT987h . HC 3 N v 7 
= 1 (v? = 223 cm" 1 ; ICernicharo et all |2000T). 1-C 3 H v 4 = 1 
(V4 = 28 cm" 1 ; Cer nicharo et al. 2000) and SiC2 v 3 = 1 (v 3 = 
160 cm-'i lGensheimer & Snvderj(l997t lCe"rnicharo et al.1 l2000). 
Radiative transfer calculations have also evidenced the impor- 
tance of infrared pumping in exciting molecules such as HC5N 
dDeguchi & Uvemurifll984l) and H 2 (lAgundez & Cernicharol 
2006) in IRC +10216. For this process to be efficient, the 
molecules must have vibrational modes active in the infrared, 
sufficiently strong and with a frequency at which the central ob- 
ject emits a high flux. The spectrum of IRC +10216 as seen by 
ISO peaks at a wavelength of ~10 fim (Cernicharo et al., 1999) 
and the flux is very large within a wide wavelength range (see 
below). In fact, the excited vibrational states of C4H, HC3N, 1- 
C3H and SiC2 detected have wavelengths between 45 and 357 
yizm. 

In the case of CH2CHCN, there exists a large number of vi- 
brational modes, the strongest of which are the bending modes 
V12 and vi 3 at 972 cm" 1 (10.3 /mi) and 954 cm" 1 (10.5 jum) re- 
spectively dCerceau et alJll985l iKhlifi et al.lfl999l) . For CH 2 CN 
the only available experim ental information conc erns the vg 
mode at 664 cm" 1 (15.1 ^m; ISumivoshi et al.|[l996l) . 

In order to give a quantitative estimate of how important is 
the infrared pumping compared to excitation by collisions, we 
may evaluate the rates of both processes. The excitation by in- 
frared pumping operates through absorption of an infrared pho- 
ton and promotion of a molecule from a given rotational level in 
the ground vibrational state (vq, J") to another rotational level in 
an excited vibrational state (vi, J), from which it decays spon- 
taneously to a different rotational level of the ground vibrational 
state (vo, J') through either a single transition or a radiative cas- 
cade process. The rate at which the latter level (vo, /') is popu- 
lated is then governed by the rate of the absorption process: vo, 
/" — > vj, J, which is given by: 



•vo,/' 



B va -> Vl 4nJ 



1 = 1 



(3) 



where B Vo ^ Vl is the Einstein coefficient for absorption and J is 
the mean intensity, averaged over all solid angles, at the fre- 
quency of the infrared band. The deexcitation from the level (vi, 
J) can be assumed as instantaneous in IRC +10216 for all rele- 
vant wavelengths. 

The rate at which the rotational level (vo, J') is populated by 
collisions with H2 molecules is given by: 



C. 



■v J' 



yuMHi) [=] s 



(4) 
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where j\ u is the coefficient for rotational excitation, within the 
ground vibrational state, by collisions with H2 and h(H_2) is the 
H2 volume density. 

To estimate the infrared flux in the outer envelope 
we have used the radiative transfer model described in 
Fonf rfa et all (120071) and the ISO observations of IRC +10216 
dCernicharoetal.L ll999). The best fit to the observed continuum 
between 7 and 27 // is shown in the lower panel of Fig. [4] The SiC 
band at 1 1 .3 fim is correctly reproduced. At long wavelengths the 
model seems to underestimate the absolute observed flux, which 
is nevertheless affected by the uncertainties in th e multiplicative 
calibr ation factors applied in the data reduction (Swinvard et all 
[19981) . 

The model has been used to calculate the radiation field at 
a distance of r = 5 x 10 16 cm (20") from the star (see upper 
panel in Fig. |4j». The calculated spectral energy distribution can 
be approximated by a blackbody with a radius = 1 3 R* and a 
temperature T^ = 550 K (see upper panel in Fig. @), so that we 
may use the Planck function for such a blackbody to evaluate 
4-nJ and see whether the rate for infrared pumping exceeds or 
not the collisional excitation rate. The radiative over collisional 
"excess" (/?_, VOj /'/C_> VOi y') may be evaluated as: 

R->v ,J' _ (gv 1 j/gv l >J") A v,^v a ( R bb/f) 2 

C^ Vo j,~ 4[«<*>W«W - \]y lu n(H 2 ) 

where (g Vu j/gv ,j") is the ratio between the statistical weights of 
the (vi,J) and (vq,J") states, v/« and A Vl _>„ are the frequency 
and Einstein coefficient of spontaneous emission of the vibra- 
tional transition, and h and k are the Planck and Boltzmann con- 
stants. Adopting, to a first approximation, (g Vu j/g Vo ,j") equal to 
unity and taking into account the variation of the H2 volume den- 
sity with radiu s, njfy) fern" 3 ! = 3.1 X 10 37 / r 2 with r expressed 
in cm (see e.g. lAgundez & Cernic haro 2006), we arrive at: 

R^-oj> = 5.8 x lQ-X^o^') 

C^ Vo J, [ e (26.2/^m)) _ l] r& ( COT 3 s -l) ( > 

The radiative over collisional excess is independent of radius, 
since both the radiation field and the H2 volume density decrease 
with the square of the radius, and only depends on the wave- 
length and strength of the vibrational band and on the collision 
coefficient. Just as an example, if we take typical values j\ u = 
10~ 10 cm 3 s _1 andA Vl _> Vo = 1 s _1 (adequate for vibrational modes 
sufficiently strong), we obtain that the rate of infrared pumping 
exceeds the collisional excitation rate for vibrational modes with 
wavelengths Air larger than ~ 6 /mi. Bending modes have usu- 
ally wavelengths larger than this lower limit, thus if the mode is 
active in the infrared and sufficiently strong, the rotational levels 
of that molecule will be mostly excited through infrared pump- 
ing in the circumstellar envelope of IRC +10216. For example, 
the two bending modes of CH 2 CHCN v 12 (972 cm -1 ; 10.3 fim) 
and V13 (954 cm , 10.5 fim) have Einstein coefficients for spon- 
taneous emission of 2.9 s _1 and 2.8 s _1 respectively The col- 
lision coefficients y u \ for CH2CHCN-H2 are not known, but ac- 
cording to equation (6), values as large as 5 x 10~ 9 cm 3 s _1 are 
required for making collisions to dominate over infrared pump- 
ing. Thus, this molecule will be preferably excited by infrared 
pumping through these two bending modes in IRC +10216. 

6. Discussion 

The detection in IRC +10216 of the partially saturated organic 
molecules reported in this Paper invites for a comparison to 

3 See e.g. http://www.lesia.obspm.fr~crovisier/basemole/ 



F„ (Jy) 




Fig. 4. Spectral energy distribution of IRC +10216 at infrared 
wavelengths. The lower panel shows the fit to the continuum 
observed by ISO. Most of the discrepancies between model and 
observations in the 7-20 fim range are likely due to molecular 
infrared bands, which are not considered in the model. In the top 
panel the radiation field (4nJ) at a point 5 10 16 cm away from 
the star, as given by the model, is compared to that expected with 
a central blackbody of radius 13 R* and temperature of 550 K. 

TMC-1, the other Galactic source that displays, together with 
IRC +10216, the largest wealth in unsaturated carbon chain 
molecules. For that purpose we have represented in Fig. [5] the 
abundances of various carbon chains in these two sources as a 
function of the number of heavy atoms. 

Chemistry proceeds in a different way in these two objects. 
In IRC +10216 interstellar UV photons drive the chemistry by 
photodissociating the closed shell molecules flowing out from 
the inner envelope. In TMC-1 the gas is well shielded against 
the interstellar UV field and the built-up of molecules is driven 
by ionization by cosmic rays. This introduces a sharp differ- 
ence in the time scales of the processes being at work. In IRC 
+ 10216 the gas travels throughout the CSE during a few thou- 
sands of years and the chemical processes take place in shorter 
time scales. In TMC-1 the built-up of molecules has time scales 
of the order of millions of years and the depletion from the gas 
of some species, due to condensation on grains, is a major issue 
at late times. A consequence of this is the lower level at which 
molecules are present in TMC-1 compared to IRC +10216. As 
an example we see in Figj5]that the most abundant carbon chains 
in both sources are those with a highly unsaturated character, 
cyanopolyynes and polyyne radicals, the abundances of which 
are larger in IRC +10216 than in TMC-1 by more than one or- 
der of magnitude. 

Another important difference between these two sources 
could lie in the C/O ratio. In IRC +10216, the processes of 
dredge-up related to the AGB phase have resulted in a pho- 
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Fig. 5. Abundances of carbon chain molecules in IRC +10216 and TMC-1. The diagram is an extension of that published by 
ICernicharo et al.l {1987). The fractional abundances relati ve to H2 are computed from the molecular column densities and the total 
H 2 column density. We use N(H 2 ) = 10 22 cm 2 in TMC-1 dCernicharo & Guelinl[l987h and N(H 2 ) = 2x 10 21 cm 2 in IRC +10216. 
The latter value corresponds to the total H 2 column density contained in an outer shell extending from 2 x 10 15 cm to 7 x 10 16 cm, 
where all the molecules considered in the diagram are most probably present. The choice of these boundaries is not critical for the 
values of the molecular abundances. Only if any of the molecules considered does not form in the outer envelope but in the inner 
regions of the CSE, then its abundance would be substantially overestimated. 



tospheric C/O ratio higher than 1 . This controls the chemistry 
in the outer envelope: almost all the oxygen keeps locked into 
CO during most of the expansion while the carbon in excess 
can participate in a rich carbon chemistry. In TMC-1 it is not 
clear whether the C/O ratio is higher or lower than 1 . The large 
number of carbon chains found argues in favor of a C/O ratio 
> 1, which is unusual for a dense cloud, and some chemical 
models have use d it because a bet ter agreement with observa- 
tions is achieved dSmith et al.ll2004l) . On the other hand, oxygen 
is cosmically more abundant than carbon and T MC-1 contains 
a substantial number of O-bearing molecules ( Ohis hi & Kaiful 
1998), which points toward a C/O ratio < 1. If the latter is true, 
then the carbon chemistry would be limited due to the presence 
of atomic oxygen which t ends to destroy the chemical com- 
plexity ( Herbs tet al -L 1 1 9941) . something that does not happen in 
IRC +10216 due to th e unavailability of free atomic oxygen 
dMillar & Herbstl [1994). This could also explain the systemat- 
ically lower abundances of carbon chains in TMC-1 compared 
to IRC +10216. 

If we focus on partially saturated species, such as 
methylpolyynes, methylcyanopolyynes, CH 2 CN or CH 2 CHCN, 
then we note that in IRC +10216 they are present at a lower level 
than the highly unsaturated species, while in TMC-1 partially 
saturated and highly unsaturated molecules have abundances in 
many cases comparable; e.g. compare the ratio of the abun- 
dances of CH 2 CN and HC3N in both sources. Thus, chemistry 
seems to favor the presence of species with a higher degree of 
saturation in dark clouds than in C-rich circumstellar envelopes. 
The recent dete ction of the highly satur ated hydrocarbon propy- 
lene in TMC-1 (Mar celino et all 20071) supports this point. 

In IRC +10216 the UV field which drives the chemistry 
is mostly photodissociating but not ionizing, so that the dom- 
inant reactions are neutral-neutral rather than ion-molecule. In 
TMC-1, the chemistry is triggered by cosmic rays ionization, 
and ion-molecule reactions do greatly participate in the built- 
up of molecules. Note that chemical models have traditionally 
explained the formation of molecules in interstellar clouds by 



means of ion-molecule reactions. However, since it was discov- 
ered that some neutral-neutral reactions are very rapid at low 
temperatures, it is now thought that these latter reactions domi- 
nate the formati on of highly unsat urated carbon chains such as 
cyanopolyynes (Herbst et al.L 1 1 9941) . That is, in general terms 
we may state that "current gas phase chemical models form 
highly unsaturated carbon chains mostly by neutral-neutral re- 
actions while partially saturated molecules are generally formed 
by ion-molecule reactions". Therefore, ion-molecule reactions 
could be, ultimately, the responsible of the higher degree of sat- 
uration observed in TMC-1 compared to IRC +10216. Grain sur- 
face reactions, which are very efficient in producing highly satu- 
rated molecules, could also play a role in the peculiar chemistry 
of TMC-1. However, it is not clear how these mantle species 
would be desorbed to the gas phase at the low temperatures pre- 
vailing in these regions. 



7. Conclusions 

In summary, we have shown that apart from the abundant 
highly unsaturated cyanopolyynes and polyyne radicals, anal- 
ogous molecules with a higher degree of saturation are also 
present in IRC +10216. In this paper we have presented the 
detection of CH 2 CHCN, CH 2 CN, CH3CCH and also of thio- 
formaldehyde. The relatively high rotational temperatures, 40- 
50 K, derived for CH 2 CHCN and CH 2 CN suggest that these 
species are excited in the circumstellar envelope through radia- 
tive pumping to excited vibrational states. 

The formation of partially saturated organic molecules in 
IRC +10216 resembles that occurring in cold dense clouds and 
stresses the similarity between the chemistry in these two types 
of sources. However, unlike in TMC-1, their abundances are 
much lower than those of highly unsaturated molecules, like 
cyanopolyynes and polyyne radicals, which reflects the differ- 
ences in the chemical processes at work in dark clouds and C- 
rich circumstellar envelopes. 
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